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Abstract 

A role of different components in the wave function of the weakly bound 
light nuclei states was studied within the framework of the cluster model, 
taking into account of orbitals "polarization". It was shown that a limited 
number of structures associated with the different modes of nucleon motion 
can be of great importance for such systems. Examples of simple and quite 
flexible trial wave functions are given for the nuclei ^Be, ®He. Expressions for 
the microscopic wave functions of these nuclei were found and used for the 
calculation of basic nuclear characteristics, using well known central-exchange 
nucleon-nucleon potentials. 
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I. INTRODUCTION. 



Theoretical consideration of strongly interacting many-particle systems, such as the 
atomic nuclei, has a number of difficulties, one of which is constructing of trial wave function. 
The wave function used should be characterized by proper set of the definite quantum num- 
bers and corresponding symmetries, taking into account the most important physical modes 
of nucleon motion. It is usually assumed that structure of the constructed wave function is 
well known in the limitsd cases, i.e. in the case where a distance between fragments is very 
large or when we have a compact one-centre system. Hence, the main attention is paid to 
determination of correct interpolation expression for the wave function in the intermediate 
region. 

Many authors have worked at simultaneous description of both "compact" and "friable" 
structures (see, for example, Ref. [0,0). Several decades ago due to activity of Yu.A.Simonov, 
great attention was paid to the hyperspherical function ||^-^ method (the K-harmonic 
method) which was considered as a general solution method for the nuclear many-particle 
problem. But it was shown that its practical application is limited effectively to compact 
systems; moreover, the results obtained hardly differ from the results calculated using 
the simplest functions of the oscillator shell model with one variation parameter. 

Recently quite effective stochastic variational methods [§|-|T^] for semirealistic nucleon- 
nucleon (NN) potentials using a number of multicluster model basic functions have been 
developed for the case of the lightest {A < 10) nuclei. They allow to find almost exact 
numerical solutions of many particle problems for some variants of NN-forces which are most 
frequently used in calculations of the light nuclei structure. Unfortunately, only numerical 
estimations of energies and radii are usually obtained due to complication of the calculations. 
Thus, results of such microscopic computations rarely can be used for estimation of other 
nuclear characteristics. Therefore, it is desirable to have many-particle wave function in 
evident form for real efficiency of the method. Moreover, if we have the wave function in the 
explicit form, we will be able to use it in calculations of different nuclear processes within 
the nuclear system. 

To construct simple and quite fiexible microscopic trial wave function of the nucleus 
which correctly describes both compact and weakly bound states, it is of prime importance 
at studing of exotic nuclei with anomalous ratio N/Z. 

It is known that deformation effects are significant in nuclei with A > A. But concrete 
calculations are usually performed on the basis of spherical symmetry for cluster wave func- 
tions and detailed function for cluster-cluster motion. Recently, we have tried ||I3| to take 
into consideration the deformation effects in the framework of the many-particle oscillator 
shell model by using the variation principle (without taking into account the cluster modes). 
In the present paper we describe microscopic model which employs orbital "polarization" as 
well as cluster degrees of the nucleon motion and exact angular projection. The model can 
be used for description of both compact and weakly bound states in the light and medium 
light atomic nuclei. 

Below we briefiy outline the construction of the trial wave function (section |I|) and give 
all matrix elements (m.e.) of the basic physical operators needed for our computation (sec- 
tion pn] ). The variational calculation results for the nuclei ^Be, ^He are given in section 0. 
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II. MANY-PARTICLE HAMILTONIAN AND TRIAL WAVE FUNCTIONS. 



Let us split considered A- nucleoli system in two subsystems {A = Ai + A2) in such 
way that nucleon-nucleon interaction inside of each subsystem is a priori stronger then that 
between the two nucleons in the different fragments. We introduce the two fragments in 
order to simplify the account. In case if there are more fragments, the generalization can 
be written in a similar way. Note that nucleon motion mode defined by interaction between 
subsystems will be the "softest". Analysis of the experimental data for the considered 
nucleus allows to choose correct splitting of A-nucleon system into fragments. As a rule, 
the most interesting "physics" at low energies is caused by the soft mode which should be 
the most detailed in the framework of the chosen model. For example, description of the 
A-nucleus weak-bound state will be sufficient if we were able to obtain the basic nuclear 
characteristics as well as the decay A = A1 + A2 threshold on the assumption that separated 
fragments Ai and A2 are known compact nuclei. 

In general, every microscopic model is associated with its own effective NN-potential. 
The most known effective NN-potentials have been constructed in the framework of the 
models based on the shell model wave function with one variation parameter. It is clear 
that using such NN-potentials in the more advanced models can be uneffective. However, 
this usage can be quite correct if NN-potential was built on the basis of the compact nuclei 
properties (e.g., ones of the magic nuclei '^He, ^^O or "^^Ca). The results of calculations [0] for 



the magic nuclei ^He, ^^O are rather similar for both simple and advanced model, therefore 
is valid to use mentioned above NN-potentials for the "polarized" orbitals model. Let us 
further consider the standard microscopic nuclear model based on the many-nucleon 
Hamiltonian with two-particle central exchange interaction 

^iji^) = ^25+l,2T+l(?^) ■ PsT{ij) , r = \fi~fj\ , (1) 

s',r=o,i 

where the radius dependence of the potential components can be represented by Gaussian 
superposition, i.e. 



,2T+1 ■ ^-^P I 



(2) 



projection operator P in eq.([l|) cut out the state of nucleon pair (ij) with defined spin S 
and isospin T. 

When fragments are separated enough, the Hamiltonian Ha should turn into the sum 
of separated subsystem Hamiltonians plus Coulomb repulsion between them due to the 
short-ranged nature of the nuclear forces 



\Rx-R2\ A 



i=\ 



i.e., the limiting cases are known. 

The behaviour of the system in the intermediate range where altitude of vector S = R\ — 
R2 significantly differs from zero, but is not large enough to neglect the nuclear interaction 
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between particles in the different clusters, is open question. Hence, it is necessary to define 
dependency of the center-of-mass kinetic energy operator Tcm. on the stage of separation of 
the A-nucleus into fragments. We simulate Tcm. in this range by expression 



(4) 



i=l 



^2 r M 



2Am 



Kj=Ai+l 



Am 



as) ■ 



E E ■ 



i=l j=Ai+l 



In our calculations the value of the parameter as {[0 < as < 1]) in eq.@ was defined by 
the degree of importance of the Pauli principle in the A-nucleon system. Thus, = if 
antisymmetrization acts on the variables of all A nucleons and = 1 in the other limiting 
case when the Pauli principle acts only on the nucleons inside the clusters. 

A total trial wave function is represented by superposition of a small number of 
determinant functions 



(5) 



where \&,y-components (normalized, but generally not orthogonal) are well known Slater 
determinants 



det 



fu 



1 



det 



<€W> 



(6) 



filled by one-particle orbitals which have unificated form: 



(7) 



Nin;R/a) ^ ^ . F ^-n, -n + ^O; ^) , 
where Xcrr is the spin-isospin function and F{a, (3; 7; z) is the Gauss hyperspherical function 

El- 

Every component \E'j, in eq.(^ represents some characteristic structure or configuration 
of the considered system and is associated to some type of nuclear motion. We hope that the 
wave function constructed in such way will be flexible enough at the reasonable number 
of the variation parameters {Si, Ri}. 
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III. OVERLAP INTEGRALS. 



In order to perform the variational calculations it is necessary to obtain the matrix 
elements of some physical operators. Needed overlap integrals can be expressed in terms 
of "partial" m.e. on the basis of the function set (|^) using known determinant function 
technique [jl6 . 



At the beginning, we consider overlap integral of the functions \E'j, with unity, i.e. 



= fufi, ■ det [{B[^^^^)ij\ 



< ^^1^^ >= /^^-det 
where the matrix is defined by 



Evident expression of the partial m.e. 

oo 

S'nma «2-R2) = J i^riA^] 0-1, Rlji^nii^'^ 0.2, R2)dx 



(8) 



(9) 



(10) 



is given in the Appendix (see eg. (|A1|)). 

The general expression for one-particle operator matrix elements 



< 



^ i=l *,i=i 



allows to obtain: 

1) m.e. for nuclear matter distribution p{f) 



1 



< 



(11) 



^AY.6{r-r..)\^,>=Y,i^i!ir-a,,R.)4^^{f-a,^ ■ (12) 

1=1 i,j=l 



2) m.e. for kinetic energy operator 



< 
where 



2=1 l,j = l 



, (13) 



< i^nMi'n, >=< i^'nli^] Si, Ri)\ - ^^l4';^^(^; Rj) > 



(14) 
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The analogous expressions can be written for < ipnMyl'ipnj > (substitution x ^ y; a ^ b) 
and for < ipni\tz\'ipnj > (substitution x ^ z; a ^ c). The matrix elements tn'n{ 

cl' , R' ] Q, -R) 

employed in eq.(|l^ are given in the Appendix (see eq.([A2|)). In eqs.(pn])-([T3|) and below 
the inverse matrix 



is used. 



Further, let us use the general formula [jT6[ for two-particle operators Yli j=i 



X 



k',l',k,l=l 



(15) 



In the case of the central exchange NN-potential (|lD, (0) we can write for even-even nuclei 
after summing up on the spin-isospin variables 



«J=1 



Z Z ^pot 

E EE<4:l(i)4:!(2)iexp 

i',i=lj',j=l 1=1 



14^1(1)4^1(2) >x 



^""^^^^ • ("^'-0. (^lii),, -^^^(0 ■ (i^r.l,)^, {b^. 

In - np 



+ 



Z N Ipot 

+ E EE<4:lw<l(2)i-p 

i',i=ij',j=i 1=1 



i4':ki)4'(2) > X 



+ 



(16) 



14^1(1)4^1(2) >x 



iV N Ipot 

+ E EE<<l(l)<l(2)|exp 
i',i=ii'j=i 1=1 



In - f2\ 



x{^'^(0 ■ (^[i]),, (^[i])^,, -^e.(/) • 



where 



3V|) + V/P 



-^ex(0 



3V^) - 



(17) 
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Y^^l^ ^ + 34'^ + VS + Vn ^ y^l) = 3ViP - 3^'^ - vS + vl? ^ 

Due to the Gaussian dependence of the considered potential the partial matrix elements 
between the basis functions (|^) are factorized on three terms, i.e. 

< V-n, (l)^n,, (2) I exp (2) >= 

= ji2{ni,ai'Ri,,n'^,aj'R^r, /i; n'laiRl, rf^ajR^j) x (18) 

X ji2(^f'&i'-Rf/, n^jihjiR^, fi] nfbiR^, n^^hjR^)ji2{nl,ceRl,, rij^CfRj,; fi; n^dR-, UjCjRj) . 
The evident expression for integral 

ji2{n'a'iR[, m'a2-R2j /^i naiRi, ma2-R2) = 

CO oo 

= y y i^n'ixi, a[R[)iJm'{x2; aai^a) exp |- ^^^^^ | tpnixi, aiRi)tprnix2; a2R2)dxidx2 

— oo — oo 

(19) 

is given in the Appendix (see eq.(|A4|)). 

In the context of described above calculations it is suitable to use the relation 

1 



, , -p(-;^J-(r37^' /^o^-yr^. (20) 







for evaluation of the Coulomb repulsion between protons. Thus, 



< 



' ^ «>J = 1 2'>j' = l «>J = 1 Q 

X < fe,,(l)fe.,(2)|exp|- ''''~ '''' I |fc,(l)fe^(2) > X 

(^lii),, (^lii),, - (^lii)., (^lii),,.} • 

i.e., in general, only integration over r remains in the equation (|21|) for the Coulomb energy 
(all two-particle integrals for Gaussian potential are expressed in terms of regular functions). 
The expression for center-of-mass kinetic energy has the form 

< *.|T,„.|^^ >= -— — < ^,1 5^ V^l^^ > + < ^ ViV.I^M > > (22) 



i=l i¥'j=l 
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where the first term in the brackets is described by eqs.([T3D, (0); the other can be evaluated 
using general eq.(|l5D for two-particle operators, i.e. 



1 ^ 
— — < V ViVjl^^ >= 



^ ^ 1/ 1 ^ /J 



A A 



5^ 5^ < 0lrWl0r^ X 0?Wl0f' > X (23) 



the explicit expressions for m.e. < 0i|V|02 > are given in the Appendix (see eq.([A^)). 

Next, we briefly consider projection of the determinant function on the state with defined 
angular momentum (the rotation effect in phenomenological models). It is obvious that 
considered above determinant functions "^i, are superpositions of the states \JM >. In 
order to obtain the state with defined angular momentum J and its projection M, we 
should define the projector Pi^j^ to operate on the total wave function \E'^. It is suitable to 
write the projector -P//^ in the Hill- Wheeler's integral form 

PifK = / dQDi/AQ) . R{Q) , (24) 

i.e. overlap integrals 



(25) 



should be evaluated to obtain eigenvalue of some physical operator O. Transformation of 
the spatial and spin-isospin functions at rotation of the reference frames on Euler angles is 
performed by usual way [[17|. A matrix B can be defined analogously to eq.(P): 

{B[.,,])^^^ J ^P^:;!ir;a,,R,)Rin)4^^if;a,,R,)df^ 




j / 7 \'"3 



-^-if) (0 if) (^) (f) ^ 



^ '^''P 1 ^ 



x-R-f {y-R\f {z-mf {i-R^y {y-R',? {~z-R^ 



z\2 



2di 262 2cf 2a] 2h] 2c] 

X = Xi, y = X2, z = x^ are the Cartesian coordinates in the rotated frame of reference that 



can be written using the rotation matrix || (f2)||, i.e. 



2'm 

m'=l 



^ Rm'mi^) ■ Xm' , m = 1, 2, 3. (27) 
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Integration in eq.(]27|) can be performed completely by transforming quadratic form from 
the exponent to diagonal one. 

Note that two-particle matrix elements containing four different orbitals can also be 
evaluated this way. Thus, in the case of m.e. 

< i^fiA^i; tti, Ri)i)fi,{r2] ttj, Rj)\exp | — 1 i?((])|^^,,(fi; a^/, i^j/)?/'^^, (fs; a^/, Rj>) > 

(28) 

the quadratic form is constructed using 6-dimension vector x = x{xi,X2, ■ ■ ■ ,xq) , where 

n = fi(xi,X2,X3), r2 = f2{x4,X5,Xe). 

Hence, all integrations over particle coordinates can be performed completely. Unfortu- 
nately, resulting expressions in the case of the angular projection are too complicated to be 
written here. Integration over Euler angles ^2, ^'3) have been carried out numerically, 
although it is possible that further investigations will allow to do it exactly. 



IV. NUMERICAL RESULTS. 

There are few extremely significant systems in nuclear physics. One of them is nucleus 
^Be which have unique properties and is critical for the theory. Thus, ^Be is 8-nucleon 
weak-bound system whose ground state is very narrow resonance near the decay threshold. 
Therefore, calculation of the ^Be basic spectroscopic characteristics will be serious test for 
the considered microscopic model. 

Due to the known symmetry of the nucleus ^Be we choose its wave function in form 

^^(1, 2, . . . , 8) = ■ ^1 + W72 ■ (^2 + ^3) , (29) 

where "^i, is Slater determinant (P) filled with simple orbitals: 
a) u = 1 

[^=1], f ^po{xj■,dl,0)^|Jo{yj■,al,0)^pQ{zj■,Cl,0), i = 1,2,3,4; , . 

\^o(a;,;a2,0)7/>o(|/,;fi2,0)^i(2,;52,0), z = 5,6,7,8; ^ ^ 



b) z/ = 2,3 



^[^^=2]/ -N r V^o(a;j;ai,0)V^o(l/i;ai,0)^o(%;ci,-s ■/?), i = 1,2,3,4; 

\i>o{xj;a2,0)ipo{yj;a2,0)'il>iizj;c2,R), i = 5,6, 7,8; 



■ n,; 



^[^^=3]/ -x ^ f V^o(a;j;ai,0)V^o(%;ai,0)V^o(%;ci,s ■ /?), i = 1,2,3,4; 

\'ilJoixj;a2,0)^oiyj]a2,0)^iizj;c2,-R), i = 5,6, 7,8; 

The values {Si, ci, 02, 62; ai, Ci, 02, C2, R, s} in these expressions are variation parameters 
defined by minimization of the total nuclear energy functional. 

Let us consider our system at different values of the parameter R in order to examine the 
system dependence on the distance between the fragments. The variation parameters were 
optimized for every value of R; when R> Rq the orbitals in the determinant functions '^i,=2 
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and \l/i,=3 are continuously transformed in order to obtain wave functions of two independent 
a-particles at the limit R oo, i.e. 



[u=2] 



i)o{xj\ ai, 0)?/;o(%-; «i, ^)'4'o{zj\ Ci, -s ■ i?), i 

-a{C) ■ ipoixj] ai, 0)ipo{yj] ai, 0)ipo{zj; ci,s- R) + 

• i^oixj] a2, 0)i>o{yj] a2, 0)?/'i(^j; C2, i?), i 

?/'o(a;j; ai, 0)V'o(z/i; «!, 0)^/'o(2:i; ci, s ■ i?), i 

«(0 ■ V'ola^j-; oi, 0)ipo{yj; ai, 0)^o(%; Ci, -s ■ R) + 

■ i'oixj-, 02, 0)ipoiyj; ^2, 0)V'i(2:j; C2, -R), i 



1,2,3,4; 



5,6,7,8; 



1,2,3,4; 



5,6,7,8; 



(33) 



(34) 



«(0 



The functions a{^), /?(^) are defined by 



i + 2ev^r^ 



/3(0 



1 



i + 2ev^r^ 



0<^<1. 
(35) 

Choice of Rq and r^o in eq.(^6|) accords to choice of interval [Rq, (?7o + l)-Ro] where one of 
the a-clusters "loses" its excitation. Next we define the antisymmetrization of the wave 
function, if R < Rq then the Pauli principle acts on all A nucleons in the nucleus. In the 
range [Rq, (r^o + l)Ro] the antisymmetrization between particles in the different clusters is 
gradually "switched off". Moreover, both "losing" of the excitation and "switching off" of 
the Pauli principle action are simulated by the same functions a{^), /?(0- 

The total wave function of the system in the considered range [Rq, (?7o + l)-Ro] is presented 

by 



vl>^(l, 2, . . . , 8) ~ /?(0 ■ + [1 + a(0] ■ ^2 ■ (^2 + ^ 



(36) 



where w^, W2 are weights of the according determinant functions aX R = Rq. The parameter 
as in the microscopic Hamiltonian (see eq.(^)) can be defined by 



as = as{R) 




R < Rq , 

Ro<R<ivo + l)Ro , (37) 
R>{Vo + l)Ro . 

In order to evaluate the overlap integrals corresponding to some physical quantity it is 
necessary to calculate the matrix B\^^^^] and the inverse matrix B^^^y In the case of nucleus 
^Be these matrices are in general 




















































































































































































det 



^4 

'50 



(38) 
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Co — OiufiPufj, ^ufi^ufi ) 



Pu^./Co -7.m/Co 

M(o -7.^/Co 

-7.^/Co 

/3,^/Co -7.^/Co .^qi 

-Su^./Co «,^/Co • ^ ^ 

-5.^/Co a.JCo 

-5.^/Co a^JCo 

-5.^/Co a.jCo 

Evident expressions for oty^^ Puimi ^v^i ^u/i can be obtained, by using the definitions ( |30| ) ( |34[ ) 
(see eqs. (|A6|) in the Appendix). 

Concrete calculations on the basis of the above considered equations have been performed 
for the case of some well known variants of the central exchange NN-potentials. Thus, the 
results for Volkov's NN-potential [|l^ (the 1st variant case) are accurate. Some of the results 
are given in Table |I[ There are energy Eq+ of the ^Be ground 0+-state, energy Ucoui of the 
Coulomb repulsion between protons, r.m.s. matter radii < >^/^ as well as optimized values 
of the variation parameters {at, bi, Ci, ai, bi, Ci, R, s} in the determinant functions and its 
non- normalized weights {wi}. In the case of Nl the wave function (|30|) contains 10 nonlinear 
variation parameters and describes both orbitals "polarization" and cluster structure of the 
nucleus ^Be. In the case of N2, the superposition of two determinant functions 

= (1, 2, . . . , A) = ■ ^i({f,}; Si = 6i, gi; ^2 = 62, ^2)+ (40) 

+W4 ■ ^4({n}; ai = bi, ci; 0.2 = C2). 

taking into account only orbitals "polarization" has been used. Column N3 corresponds to 
the previous case with neglecting of the deformation, i.e. 

^'a = (1, 2, . . . , A) = ■ *i({fi}; S) + ■ *4({r=;}; ~a). (41) 

The values of column N4 are calculated on the basis of one determinant function with 
different degree of "polarization" for s- and p-orbitals (4 variation parameters: fii = 61, 
Ci, 02 = 62; C2). In the case of N5 we have used one determinant function with the same 
deformation for s- and p-orbitals (2 parameters: a = b, c). The simplest approach is 
microscopic SU(3)-model with one variation parameter a = b = c that corresponds to 
column N6. The data in column N7 have been calculated using the most complicated wave 
function It differs from one of the case Nl by additional shell- type determinant function 

'^u=4: filled by one-particle states containing independent variation parameters Si = 61, ci, 
02 = b2, C2- This approach employs 14 nonlinear variation parameters. And, finally, data 
in the column N8 have been evaluated using the previous function "^a, but without angular 
projection. 

First we consider numerical results obtained with one-center functions \I'^ (models N2- 
N6). The data in Table | show that the most binding energy growth is caused by taking into 
account the oscillator field deformation (compare models N6 and N5). Independent orbital 
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"polarization" is also significant (compare models N5 and N4). If the nonlinear variation 
parameters {at, bi, Cj} are the same (i.e. deformation is neglected) than expanding of the 
basis from 1 to 2 determinant functions \E'jy does not practically change the nuclear energy E 
(comp. models N3 and N6). Moreover, our estimations show that the employment of two, 
three and more determinant functions with the same SU(3)-symmetry does not allow to 
reach significant advance even if all considered s- and p-orbital deformations are taken into 
account (comp. models N2 and N4). Hence, the deformation effects can be considered using 
only one one-center Slater determinant. More effective is expansion of the variation basis by 
cluster wave function. It enhances the '^Be nuclear energy by approximately 3 MeV (comp. 
models Nl, N4, N2). 

Fig. ID shows ^Be mass distribution p(f^ evaluated by using the wave function (|30|). The 
optimum variation parameters correspond to model Nl (see Table |I|). The calculations are 
carried out in the nuclear center-of-mass reference frames. Surface {x, y = 0, z) is displayed 
in Fig. |I| due to axial symmetry of the distribution p{r). We can see from Fig. |I]that distance 
between density maximums (centers of mass for two a-particles) is D ~ 4.2 fm. Density 
of junction between a-particles is approximately 0.7pmax- The distribution p(r) obtained in 
the framework of the other considered models is rather similar to Fig. |l|, although distance 
D and width of junction can be essentially different. 

Evaluated parameter R dependence of the ^Be total energy functional /e^R, h, Cj}) is 
plotted in Fig. 0. The shown dependence have been calculated under condition that for every 
value R functional Je is optimized to the other variation parameters. Such computations 
can be useful for estimations of the nucleus-nucleus potentials. Note that in a range where 
clusters overlap each other, the distance D between density maximums can essentially differ 
from its asymptotic value which is equal to Das = -R(l + s) in our case. Fig. |^ shows 
dependence D = f{R) obtained for model Nl. We can see that the distance D can slightly 
decrease at < i? < 1.5 fm and only at i? > 2 fm the R dependence of D is almost linear. 

Table |T| contains characteristics of the ^He ground 0"'"-state calculated using three vari- 
ants of NN-potentials: VI from Ref. 0, BBl from Ref. [H and G5 from Ref. 0. Com- 



putations are carried out in the framework of two models. Wave function of model I is a 
Slater determinant for 6 nucleons filled by orbitals (|3lD, and contains four variation param- 
eters: ai = bi, Ci for first four orbitals and 0,2 = &2; C2 for the last two orbitals. Wave 
function of the second model is superposition of a few Slater determinants, similarly to the 
^Be: one-center determinant (model I) and two-center determinants with centers in points 
Ri = (0, 0, ±i?) and R2 = {0,0, ^sR). This model contains 10 nonlinear variation param- 
eters whose optimum values for the considered variants of NN-forces are given in Table 

El 

It should be remembered, that a core in ^He is nucleus ^He. Its energy calculated using 
simple SU(3)-model function is equal for three variants of NN-forces: £^a(Vl) = —27.09 
MeV, ^a(BBl) = -27.37 MeV, E„(G5) = -28.29 MeV. Comparison of these values with 
^He ground state energies in Table |l| shows that two neutrons are bound only in the case 
of model II and NN-potentials VI and G5. Therefore, we can conclude that wave function 
of model II describes the real systems more accurate than one-center function of model I. 

Note that authors don't claim superexact description of the considered here nuclear 
systems which has been already studied profoundly by many authors (see, for example, 
Refs. P,p!0[]). The purpose of this study is to clarify the question about significance of 
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"compact" and "loosely bound" structures at consideration of nuclear states using simple 
trial wave functions. 

Performed analysis allows to draw up some general scheme for the construction of the 
simplest and quite flexible trial wave function of weakly bound state as well as a strategy 
for concrete computation. First, possible significance of the deformation in formation of 
considered nucleus state is roughly estimated in the framework of the one-center determinant 
function model. Next, cluster degrees of freedom are taken into account. Hierarchy of 
the nucleon motion modes can be determined using non-projected wave functions that 
significantly reduces computations. Moreover, such reduced evaluations give approximate 
optimum values for almost all variation parameters (for example, compare corresponding 
parameters of models N7 and N8 in Table |). It is obvious that this fact is especially 
significant in order to find the global minimum of the total nuclear energy functional by 
employing projected functions 2, . . . ,A). 



APPENDIX: 



The evident expression (0) allows to calculate: 
1) partial m.e. ([T0|) 



5'nin2('^l; Rll 0,2, -^2 



9o 



N (n,; f ) N (n. 



a2 



ni+n2 



(Al) 



xF 



ni+n2 ni+n2 I 1 
2 ' 2 ^ 2' ' Z| 



90 



= a/AqIo ■ exp 



(Ri - R2) 
2{al + al) J ' 

1 [Ri 



An 



2 1 2 ' 



^0 — 77 ( — ^0 H ^0 ) ; 

2 V«i ^2 



2) m.e. ([T^) used for evaluation of the kinetic energy 



tn'Tiia', R!\ a, R) =< 'ipn'i.x; a', R' 



2m dx'^ 



2a\ 



2 I 2 ' 



\^n{x\a,R) >-- 



(A2) 



+2n—N in - 1, — ] gn\n-i{a', R'; a, R) + ( ^ - 2n - l) N ( n, — ] gn',n{a\ R'; a, R)- 



-2^N (n + 1, ^ j gn',n+i{a', R'; a, R) + N (n + 2,^ ] gn>,n+2{a', R'; a,R)}; 
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d d d 

< (/'i|V|02 >=< (t)i\-Q^(t)2 > ■ex+ < <^il^<^2 > ■ey+ < h\-Q^<f>2 > -ez , (A3) 



d ^ ~ d 



Ry-hRy)gn^nXc,R'\cR^)^ (A4) 



X 



tT. TDX. „ -DX\ 

n ^T-, -. — — —gnXn^_i[a,R ,aR ) + —gn^^n^[a, R ,aR ) — 



A^(n^ + l;i?^/a) 

> and < 01 1 
substitution a ^ b, R^ ^ R^ and a ^ c, R^ ^ R 



5'n^,n^+l(C'; R^', aR^ 



equations for < 0i|^02 > and < 0i|£02 > can be obtained from ( |A3|) after obvious 



3) 2-dimension integral (|T^) determining the potential energy of NN-interaction 



ji2{naiRi,ma2R2', naiRi, ma2-R2) 



1 



(A5) 



2ai/ \ ai V2«2/ \ 02 / V2ai/ V '^1/ V2a2 



I m; — 1 X 

(32 



X 



\ fi+n / ^ \ rh+m n+n //- , m [2] 

) ("I) -(-->E^E 



n + n-/)! ^ A;!(/-2A;)! 



■X 



X 



/-2fe 



[ n + TT- + "T' + y»' — ^ 1 



i=0 



(ra + n + m + m — /) ! 
i!(f;, + r;, + m + m — / — 2z)! 



where 



X 



/i / R\ R\ \ /i^o / -^2 



- h+n+fh,+m—l—2i 



+ ^ + 



^0 



1 + — — h 
2al 2al 



(A6) 



2al 2~al 2al 2~al 



2^^ 25f 



2^ ^ 2S| 



I T)2 p2 p2 p2 

1 / it]^ iX]^ 1X2 

2 \ af a\ a\ 



+ 
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Due to definitions the B[^^^] matrix elements (see eq.(PBD) can be expressed in 

terms of functions gnm from eq.(^): 

"11 = "22 = 1, "12 = "21 = fi'oo(^i' 0; ^11' 0)5^00(51, 0; Cii, 0), 

/5ii=/522 = l, /5i2 = /521 = 5'oo(«2,0;a22, 0)5-11(02, 0;c22,0), 

711=722 = 0, 712 = 721 = fi'ool^i'O; 022, 0)501 (ci,0;c22,0), 

Sn = 622 = 0, 612 = S21 = 5oo(^2, 0; an, 0)510(02, 0; cn, 0). 
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FIGURES 



FIG. 1. Mass density of ^Be as a function of the variables z and r_L = \J -\- y^, and contours 
of constant mass density in the section of ^Be by the (x, z) plane in the intrinsic coordinate frame 
rigidly bound to the nucleus. 

FIG. 2. Total energy of the ^Be nucleus as a function of the variational parameter R. 

FIG. 3. 2Q;-distance as a function of the parameter R. 
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TABLES 



TABLE L The nucleus ^Be ground state in microscopic models, Volkov's NN-potential |16|. 
The models (N1-N8) are described in the text. 



Characteristics 


1\T 1 


M 9 


1\T 

IN 




1\T K 
IN 


IN D 


M 7 


M S 
L\ 


E 


MeV 


CO 700 


-DU.Ui ( 


-oD.oDy 


AO (ir\K 


/1 7 01 n 
-4 / .yiu 


-oD.loO 


CO Q7c; 
-Oz.y ( 


AA SQK 
-44.00D 


Ucoui, MeV 


Q 1 /1 1^ 
0.140 


o.oou 


ij.4ZZ 


0.006 


o.oDD 


0.44D 


O.i04 


Q 07Q 

o.z /y 


, 2 


^1/2 r 

>^'^, fm 


z.oz 


Z.oo 


1 Q 

ly 


z.oi 


Z.oU 


Z. 10 


z.Oz 


z.oo 




ai — bi 


l.Zl 




1 /in 


1 on 
i.zy 


1 QO 
l.OZ 


i.uU 


1 00 
l.zz 


1 07 
l.z / 


wi 


Cl 




z.oo 


1 4Q 


9 4"^ 


9 n? 


1 fin 


9 1 R 


1 SQ 

-L .Oc/ 




0,2 = 62 


1.35 


1.29 


1.49 


1.35 


1.32 


1.60 


1.29 


1.34 




C2 


1.91 


1.80 


1.49 


1.93 


2.07 


1.60 


1.79 


1.62 




Ol = 


1.49 












1.47 


1.57 




Cl 


1.59 












1.56 


1.56 




02 = 62 


1.28 












1.25 


1.30 


^3 


C2 


1.37 












1.36 


1.36 




S 


1.28 












1.29 


1.30 




R, fm 


1.76 












1.73 


1.53 




ai = 61 




1.50 


1.86 








1.50 


1.51 




Cl 




2.70 


1.86 








2.48 


2.10 




0.2 = 62 




1.53 


1.86 








1.65 


1.72 




C2 




2.28 


1.86 








2.37 


2.12 




Wl 


1.74 


1.35 


1.03 


1.00 


1.00 


1.00 


0.97 


0.50 




= W3 


1.00 












0.69 


0.21 




W4 




0.72 


0.69 








0.42 


0.34 
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TABLE II. Calculated ^He features. Ener gies E, Ucoui are given in MeV, rms radii, {oj, 6j, Cj}, 
R are given in fm. 



Poten. VI 



Poten. BBl |18] 



Poten. G5 |19| 



Characteristics 


I model 


II model 


I model 


II model 


I model 


II model 


E 


-25.296 


-27.396 


-24.451 


-26.918 


-29.636 


-30.831 




0.7139 


0.7633 


0.6820 


0.7428 


0.7300 


0.7591 


matter 

TTYIS 


2.51 




2.64 




2.47 




proton 
frryis 


1.86 




1.95 




1.78 




^neutron 
rms 


2.78 




2.92 




2.76 




di = bi 


1.31 


1.28 


1.36 


1.34 


1.34 


1.33 




2.19 


2.05 


2.32 


2.16 


1.97 


1.85 


do = bo 


1.45 


1.36 


1.55 


1.47 


1.60 


1.74 


C2 


2.14 


2.05 


2.23 


2.16 


2.00 


1.83 


Ol = 61 




1.39 




1.41 




1.39 


Cl 




1.37 




1.39 




1.41 


*i/=2,3 0,2 = 62 




1.86 




1.85 




1.46 


C2 




2.54 




2.56 




2.47 


S 




1.16 




1.13 




1.65 


R 




1.31 




1.44 




0.91 


Wl 


1 


0.9209 


1 


0.8509 


1 


0.8993 


W2 





0.6223 





0.7033 





0.4835 
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R [fin] 



9- 




R [fin] 



